Gamma-ray flares from the Crab Nebula 
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A young and energetic pulsar powers the well-known Crab Nebula. Here we 
describe two separate gamma-ray (photon energy >100 MeV) flares from this 
source detected by the Large Area Telescope on board the Fermi Gamma- 
ray Space Telescope. The first flare occurred in February 2009 and lasted 
approximately 16 days. The second flare was detected in September 2010 and 
lasted approximately 4 days. During these outbursts the gamma-ray flux from 
the nebula increased by factors of four and six, respectively. The brevity of 
the flares implies that the gamma rays were emitted via synchrotron radiation 
from PeV (10 15 eV) electrons in a region smaller than 1.4 x 10 2 pc. These are 
the highest energy particles that can be associated with a discrete astronomical 
source, and they pose challenges to particle acceleration theory. 

The Crab Nebula is the remnant of an historical supernova (SN), recorded in 1054 C.E., 
located at a distance of 2 kpc (/). The SN explosion left behind a pulsar, which continuously 
emits a wind of magnetized plasma of electron/positron pairs (henceforth referred to as elec- 
trons). This pulsar wind is expected to terminate in a standing shock where the particles may 
undergo shock acceleration (2, 3). As the electrons diffuse into the downstream medium they 
release energy through interactions with the surrounding magnetic and photon fields. This emis- 
sion is observed across all wavebands from radio up to TeV gamma-ray energies and is referred 
to as a pulsar wind nebula (PWN). The efficiency of this process is remarkable. As much as 
30% of the total energy released by the Crab pulsar is emitted by the PWN [ (4) and references 
therein]. The Crab PWN has an approximately ellipsoidal shape on the sky with a size that 
decreases with increasing photon energy. At radio frequencies it extends out to 5' (3 pc) from 
the central pulsar. At X-ray wavelengths a bright torus surrounds the pulsar; its radius is 40" 
(0.4 pc) and jets emerge perpendicular to it in both directions. 

Within the region encapsulated by the torus there are several small-scale structures. The 
inner nebula, which we define as the central 15" around the pulsar, has several small-scale 
regions of variable X-ray and optical brightness. The most prominent is an X-ray-bright inner 
ring with a radius of 10" (0.1 pc); this ring is thought to represent the termination shock of the 
PWN (5). Several knots with diameters of ~1" (0.01 pc) are detected close to the inner ring and 
the base of the jets, and bright arcs of comparable width are observed moving outwards from 
the inner ring into the torus (6, 7). 

The broad-band spectral energy distribution (SED) of the Crab Nebula is composed of two 
broad non-thermal components. A low-energy component dominates the overall output and ex- 
tends from radio to gamma-ray frequencies. This emission is thought to be from synchrotron 
radiation. This notion is confirmed in radio to X-ray frequencies with polarization measure- 
ments (8-10). The emission of this synchrotron component peaks between optical and X-ray 
frequencies, where the emission is primarily from the torus (5). The emission site of higher- 
energy photons (beyond 100 keV) cannot be resolved due to the limited angular resolution of 
telescopes observing at these frequencies. The high-energy component dominates the emission 
above ~ 400 MeV and is thought to be emitted via inverse Compton (IC) scattering, predomi- 
nantly of the synchrotron photons (11, 12). 
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The large-scale integrated emission from the Crab Nebula is expected to be steady within 
a few percent and is thus often used to cross-calibrate X-ray and gamma-ray telescopes and to 
check their stability over time (13, 14). Recently, variability in the x-ray flux from the nebula 
by ~3.5% yr _1 has been detected, setting limits on the accuracy of this practice (15). Yearly 
variations in the emission in the high-energy tail (1-150 MeV) of the synchrotron component 
has also been reported (16,17). No significant variations have been detected for the high-energy 
component of the nebula (18-20). 

The Large Area Telescope (LAT) on board the Fermi Gamma-Ray Space Telescope (Fermi) 
has continuously monitored the Crab Nebula as a part of its all-sky survey since August 2008. 
The LAT detects gamma rays from 20 MeV to >300 GeV, and this spans the transition region 
between the low and the high-energy components of the nebular spectrum. The average SED 
measured during the first 25 months of observations (Fig. 1) is well characterized by the sum 
of two spectral components, each with a power-law dependence on energy (21). The integrated 
flux of the low-energy component is (6.2 ± 0.3) x 10~ 7 cm~ 2 above 100 MeV with an photon 
index of 3.69 ± 0.11 (only statistical errors are given; see online supplements for a discussion 
of systematic errors). The high-energy component has an integral flux of (1.3 ±0.1) x 10~ 7 
cm -2 s _1 above 100 MeV with a photon index of 1.67 ± 0.04. Due to its hard energy spectrum 
the high-energy component dominates the emission above 426 ± 35 MeV. 

In order to search for flux variability of both spectral components in the LAT band, we 
grouped the flux measurements into monthly time bins. The high-energy component was found 
to be stable. The low-energy component was found to vary on these time scales (Fig. 2); 
the probability that the measured flux variations are statistical measurement fluctuations in a 
constant source is less than 10~ 5 . No significant spectral variations were detected for either 
component on monthly time scales. Flux variability was also searched for on sub-monthly time 
scales, for which the low-energy component of the nebula is significantly detected by the LAT 
only in high-flux states. The flux of the low-energy component was significantly enhanced 
compared to the average values in February 2009 and September 2010 (Fig. 2). No variations 
were found for the high-energy component. The September flare was first announced by the 
AGILE gamma-ray mission (22), which additionally reports a flare in October 2007, before 
the start of Fermi observations (23). The Fermi-LAT detected flare in February 2009 was not 
detected by AGILE as the instrument was pointing at a different part of the sky. 

The February flare had a duration of ~16 days. The average integral flux above 100 MeV 
of the low-energy component between MJD 54857.73 and 54873.73 was (23.2 ± 2.9) x 10~ 7 
cm~ 2 s -1 , corresponding to an increase by a factor 3.8 ± 0.5 compared with the average value; 
the increase is significant at >8-er level. The September flare lasted for only ~ 4 days. The 
integral flux above 100 MeV between MJD 55457.73 and 55461.73 was (33.8 ± 4.6) x 10~ 7 
cm -2 s~\ corresponding to an increase by a factor 5.5 ± 0.8 with respect to the average and a 
significance of >10 a. 

The February flare has a soft spectrum with a photon index of 4.3 ± 0.3 (Fig. 1). The 
spectral slope is compatible with the average 25-month value within two standard deviations. 
The energy spectrum for the second flare was significantly harder, with a photon index of 2.7 
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± 0.2, and was still detected above 1 GeV at a 3-a-level. The average power released in each 
of the gamma-ray flares was approximately 4 x 10 36 erg s _1 , for the case of isotropic emission. 
No significant variations in the emission of the pulsar were detected on monthly and four-day 
time scales through the period of observations. Examination of the timing residuals of the 
pulsed emission indicated no significant variations during either flares nor any significant glitch 
activity during the first 25 months of LAT observations. 

No variations in the synchrotron component between infrared and X-ray frequencies were 
seen about the average nebular flux level during the second flare (24). We analyzed data col- 
lected by the BAT instrument on board the Swift satellite (25), which continuously monitors the 
sky at photon energies of 15 — 150 keV. The mean flux measured during the first flare was (2.0 
±0.1) x 10~ 8 erg cm~ 2 s _1 , the flux during the second flare was (2.0 ±0.1) x 10~ 8 erg cm~ 2 
s _1 . Both observations are therefore within 5% of the average flux of (2.09± 0.10) x 10~ 8 erg 
cm -2 s -1 measured by BAT in this energy range (2(5), and show no correlation to the gamma- 
ray flares. The angular resolution of the BAT only allows for the measurement of the spatially 
integrated spectrum. Sub-arc-second resolution images were taken in X-rays by the Chandra 
observatory and optical by the Hubble Space Telescope a few days after the second flare. Al- 
though both images show no unusual activity compared to previous observations, both show a 
brightening 3" east of the pulsar (27). In the Chandra image this brightening is associated with 
a knot of ~1" diameter that might be associated with the inner ring or the base of the jet. Such 
a brightening might be interpreted as an afterglow at lower frequencies of the gamma-ray flare, 
but no conclusions can be drawn based on one event. 

The brief flare time scales and the requirement that the emission volume be causally con- 
nected imply that the flaring region must have been compact. If L is the diameter of the flaring 
region along the line-of-sight and t is the flare duration, then L < Dct, where the Doppler fac- 
tor D accounts for relativistic boosting effects. The Doppler factor is expected to be moderate 
within the Crab Nebula, as the typical velocities observed are smaller than 0.9 c (7). Even if the 
emission region was moving directly toward us, this yields D <4.4. For a flare duration of 4 
days this results in L < 1.4 x 10~ 2 pc, which corresponds to < 1.5 " projected on the sky. Struc- 
tures this small are found only in the inner part of the nebula, close to the termination shock, 
the base of the jet or the pulsar, suggesting that the gamma-ray emission detected in the flare 
originated from these regions. This is in agreement with expectations of relativistic magneto- 
hydrodynamic simulations, in which the gamma-ray emission of the synchrotron component 
originates close to the termination shock (28, 29). 

The extrapolation of the the LAT spectrum of low-energy component to lower frequencies 
suggest that it represents synchrotron emission (Fig. 1). The brevity of the gamma-ray flares 
strengthens this scenario: If the flare were instead produced by IC radiation or Bremsstrahlung, 
the cooling time of the emitting electrons would greatly exceed the flare duration. The cooling 
via Bremsstrahlung in particle densities <10 cm -3 (30) happens over ~10 6 years. Similarly, 
electrons cooling via IC emission of 100 MeV gamma rays on the photons of the synchrotron 
component of Crab Nebula have cooling times > 10 7 years. The average magnetic field inside 
the Crab Nebula is estimated to be ~200 /iG, as deduced from modeling of the broad-band 
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SED (12,21), and might be enhanced locally by up to an order of magnitude in the inner neb- 
ula (31). These fields imply synchrotron cooling times <15 days, comparable to the flare dura- 
tion, leaving synchrotron radiation as the only plausible process responsible for the gamma-ray 
emission during the flares. 

The detection of synchrotron photons up to energies of > 1 GeV confirms that electrons 
are accelerated to energies of >1 PeV in the Crab Nebula (32). These are the highest energy 
particles that can be associated directly with any astronomical source, and they pose special 
challenges to particle acceleration theory. Because synchrotron losses are so efficient, there 
must be a strong electric field E to compensate radiation reaction, given by: 

E/B « r L /4ooi > (1.3aS lpk /m e c 2 ) « (£ 7pk /50MeV) (1) 

where r L is the Larmor radius, £ coo \ is the radiative cooling length, a is the fine structure con- 
stant and £ 7pk is the peak synchrotron frequency at which the most energetic electrons are 
emitting (17, 33). Due to the detection of gamma-ray emission beyond 1 GeV £ 7pk can be 
conservatively estimated to be greater than 200 MeV. The electric field is unlikely to exceed the 
magnetic field; if it did, there would be a local reference frame with a pure electric field in which 
vacuum breakdown would occur quickly. We conclude that the electric field, as measured in the 
Crab frame, is close in magnitude to the magnetic field in the region where the highest energy 
synchrotron photons were emitted. This subsumes the possibility of bulk relativistic motion. 
Furthermore, the resistive force due to radiation reaction is competitive with the Lorentz force 
and the cooling length is comparable with the Larmor radius. This poses severe difficulties 
to the widely-discussed acceleration mechanism of diffusive shock acceleration (34, 35). The 
proposed acceleration due to absorption of ion cyclotron waves does not suffer from these con- 
straints (36). However, it appears to operate on time scales which are too long to accommodate 
the fast variability seen during the flares. Alternatively, the acceleration could be related directly 
to the electric field from the pulsar. 

The Crab Nebula is powered by the central neutron star which acts as a DC unipolar inductor 
and a source of an AC striped wind (2, 3). What happens to the DC and AC current flows is 
controversial. It is widely supposed that ~ 90% of the DC current returns in an outflowing 
wind that becomes particle-dominated and encounters a (mostly invisible) termination shock at 
a radius ~ 0.1 pc (37), but the wind could also remain electromagnetically dominated (33, 38). 
For the measured spin-down rate, a moment of inertia of ~ 1 x 10 45 g cm 2 and a force-free 
model of the magnetosphere, the total induced potential difference is ~ 50 PV, high enough 
to accelerate particles to the required energies. The current associated with this potential is 
~ 300 TA yielding a DC power per hemisphere of ~ 1.5 x 10 38 erg s _1 , a factor ~40 larger 
than the power released in the flares. Another interesting possibility is that particle acceleration 
takes place in the AC striped wind of the pulsar due to magnetic reconnection, although it is not 
clear if this process can accelerate particles to PeV energies on the required time scales (39, 40). 

The observations reported here have raised compelling questions on our understanding of 
particle acceleration and motivate more detailed calculations; together with the ongoing gamma- 
ray observations of the LAT and observational campaigns at X-ray and optical wavelengths they 
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might soon pinpoint the gamma-ray emission site in the Crab Nebula. 
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Figure 1: Spectral energy distribution of the Crab Nebula. Black open circles show the average 
spectrum measured by the LAT in the first 25 months of observations. Red squares show the 
energy spectrum during the flare of February 2009 (MJD 54857.73-54873.73) and blue open 
squares the spectrum in September 2010 (MJD 55457.73-55461.73). Gray squares show his- 
torical long-term average spectral data from the COMPTEL telescope with 15% systematic 
errors (41). Arrows indicate 95% confidence flux limits. 
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Figure 2: Gamma-ray flux above 100 MeV as a function of time of the synchrotron component 
of the Crab Nebula. The upper panel shows the flux in four-week intervals for the first 25 
month of observations. Data for times when the sun was within 15° of the Crab Nebula have 
been omitted. The gray band indicates the average flux measured over the entire period. The 
lower panel shows the flux as a function of time in four-day time bins during the flaring periods 
in February 2009 and September 2010. Arrows indicate 95% confidence flux limits. 
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Supplements 



Fermi-LAT Analysis 

The LAT analysis presented in this publication was performed using the Science Tools v9rl7p0 
with the P6_V3 -DIFFUSE instrument response functions (42). Fluxes were estimated by max- 
imizing the likelihood of a given source model using the unbinned gtlike tool. The model was 
composed of the Galactic and isotropic diffuse models (v02 (43)), as well as all sources within 
15° of the Crab pulsar in the first Fermi LAT catalog (1FGL (43)). Only photons within this 
same radius and with energies greater than 100 MeV were considered. All data taken under 
optimal conditions between MJD 54684.73 and 55477.73 were considered, except for periods 
when the sun was within 15° of the Crab Nebula. To avoid background contamination from the 
gamma-ray emission of the earth limb only time periods where the Crab Nebula was at a zenith 
angle <90° were considered. In order to avoid the strong foreground emission of the pulsar, 
only rotational phases between 0.52 and 0.87 were considered for the nebula analysis (21). The 
pulsar timing model used will be made available through the Fermi Science Support Center (42). 

Light curves were generated by performing the likelihood fit in each time bin independently. 
The free parameters in the 28-day likelihood fits were the fluxes and spectral indices of the 
synchrotron and IC components of the nebula and the amplitude of the Galactic and isotropic 
diffuse emission. All other parameters were held fixed to their 1FGL values. In the 4-day light 
curves, the only free parameters of the likelihood fits were the flux and spectral index of the 
synchrotron component. The IC flux and spectral index and the Galactic and isotropic diffuse 
normalizations were held fixed to their 25-month best- fit values. Spectral data points on the 
SED were measured by performing the likelihood fit in each energy bin independently, varying 
only the normalization of the best-fit model. For both light curves and energy spectra, 95% 
confidence upper limits on the flux were derived if the significance of the source detection was 
below 2a. 

Significance and position of the flares 

The flaring periods were identified by fitting a flaring component together with the average 25- 
month spectral model in weekly time bins between MJD 54684.73 and 55461.73. The average 
model included the two components of the Crab Nebula. The flaring component was modeled 
by a power-law energy spectrum and its normalization and index were the only parameters 
varied in the fits. The significance of the flaring component was <3.5 a for all weeks outside of 
the flare periods, while it was >6.7 a during the flares. Considering trials for the 102 analyzed 
weeks this corresponds to significances <2 and >6 a respectively. To measure the properties 
of the emission during the flare in more detail two flaring periods were defined from the 4 
day light curves shown in Figure 2. For the first flare the flaring period was defined as MJD 
54857.73-54873.73 (the four high flux bins around the flare). Fitting the additional source with 
the average best- fit model yields a significance of >8 a for the flaring component. For the 
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second flare the flaring period was defined as MJD 55457.73-55461.73. The significance of the 
flaring component during this time period is >10 a. 

For both flares the position of the flaring component was derived using the gtfindsrc program 
of the Science Tools. Both positions are compatible with the Crab Nebula. The best-fit position 
of the first flare is R.A.: 83.73°, Dec: 21.81° with a 68% containment radius of 0.22° (J2000), 
one standard deviation away from the Crab pulsar position of R.A.: 83.63° Dec: 22.01°. The 
best-fit position of the second flare is R.A.: 83.68°, Dec: 22.03° with a 68% containment radius 
of 0.044°, 1.1 a away from the pulsar. 

The possibility of a background blazar 

The monthly flux variability and flaring periods in principle might be due to the emission of a 
background source that cannot be resolved within the angular resolution of the LAT. However, 
this scenario is very unlikely for several reasons. First, the flux increase during both flares above 
100 MeV was >2x 10 -6 cm -2 ; there are fewer than 100 LAT sources which show variability 
of this magnitude. The probability of chance coincidence for one of them to be located behind 
the Crab Nebula within the 2 a localization error measured for the second flare, is <6 x 10~ 5 . 
Second, the spectral index during the first flare is the softest of any source yet detected by 
the LAT, with the exception of the synchrotron component of the Crab Nebula itself (21, 43), 
suggesting that the emission originate from the same source (this argument also relates the 
second flare to the Crab Nebula, as the probability that both flares are produced by two different 
sources is negligible, due to the even lower chance probability of positional coincidence in this 
case). Finally, the only source class known to produce variability on the flare time scales of 
the observed magnitude in gamma-rays are blazars. No known blazar is located near the Crab 
Nebula within the angular resolution of Fermi and X-ray observations taken two days after the 
second flare revealed no new source which might be associated with a yet-unknown blazar (44). 

Systematic Uncertainties 

Stability of flux measurements with the LAT data over time was tested using sources which 
are expected to be constant in flux. These sources are the Crab pulsar, the nearby Geminga 
pulsar, the Vela pulsar and the Galactic diffuse emission. We found the light curve of each of 
these reference sources was consistent with a constant on monthly and 4-day timescales if we 
added a 5% systematic error in quadrature with the statistical errors. We have therefore added 
a 5% systematic error in quadrature to all light curves in this publication and have taken this 
systematic error to be time-variable and normally distributed. We adopt 30% as the systematic 
uncertainty for overall normalization of integrated fluxes >100 MeV reported here. The errors 
given on LAT spectral points are statistical only. The uncertainties in the energy dependence of 
the acceptance and energy disperson of the LAT leads to a systematic error of 0. 1 on the photon 
index (21). 
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